solve d = v*cos(theta)/g*(v*sin(theta)+sqrt((v*sin(theta))^2+2*j*y)) for theta